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ESTABILIDADE AERODINÂMICA DE CHAMINÊS, 
TORRES E ANTENAS 


RESUMO 


Este artigo constitui, com ligeiras modificações, o capi- 
tulo 4 do Curso de Dinâmica de Estruturas, recentemente 
editado pelo Laboratório Nucional de Lingenharia Civil. 

Referem-se as principais caracteristicas dinâmicas das 
estruturas de chaminés, torres e mastros espiados e quanti- 
ficam-se as solicitações devidas à acção do vento sobre elas. 
Analisam-se os diversos tipos de comportamento em regime 
estático e em regime dinâmico, utilizando uma formulação 
matricial corrente, Descrevem-se os principais aspectos 
a terem conta no dimention unento de espiar, In licam-se 
finalmente processos para atenuar ou eliminar a formação 
de turbilhões de Karman, 


1 — GENERALIDADES 


A construção de chaminés metálicas e de betão 
armado de grande altura bem como de torres 
auto-resistentes ou espiadas (antenas) destinadas 
a diversos fins e em particular a telecomunicações, 
obriga a considerar de um modo especial a res- 
pectiva estabilidade aerodinâmica. 

As principais razões que concorrem para este 
facto são : 


1) As características dinâmicas deste tipo de 
estruturas : 

A frequência própria fundamental é geralmente 
inferior a 1 Hz e o amortecimento relativo 1 não 
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ARTUR RAVARA 


Engenheiro Civil. 


por 


Especialista do Serviço de Edifícios 
e Pontes, Chefe da Divisão de Dina- 
mica Aplicada, do LNEC. 


SYNOPSIS 


The present article is based on the 4th chapter of a post= 
graduate course on Structural Dynamics recently lectured 
by the author and edited by the Laboratório Nacional 
de Engenharia Crvil, 

The main dynamic characteristica of atacks, towers and 
quyed masts are indicated, The static and dynamic loads 
due to wind action are quantified. The structural behaviour 
under those actions is analysed by means of usual matrix 
formulation, Typical problems arisen by quy cables are 
briefly dealt qui, Finally practical methods to avoid 


or to aittenuate Karman vortexes are indicated, 


vai além de cerca de 3º/,, podendo indicar-se 
como valores de referência: 


— estruturas de aço soldado não revestido: 
n=0,2 à 0,5 / 


— estruturas de aço rebitado ou aparafusado: 
f=10 à 15 


— chaminés de aço soldado revestidas: 
n==1,0 a 1,5º/, 


— chaminés de betão armado e pré-esforçado: 
n=2 a 3º, 

A tendência em chaminés metálicas é genera- 

lizar o emprego de ligações por soldadura e utilizar 


383 


aços com tratamento anticorrosivo, o que dis- 
pensa revestimento interior. As duas disposições 
conjugam-se para baixar o amortecimento. 


2) Em geral são estruturas muito expostas ao 
vento pela sua grande altura e por se situarem 
em locais elevados (em especial as torres de tele- 
comunicações). 


3) O dimensionamento é frequentemente con- 
dicionado por limitações de deformabilidade : em 
chaminés as deformações podem prejudicar a 
estanqueidade de ligações a tubagens, sobretudo 
se nelas circulam produtos corrosivos ou infla- 
máveis; em torres e antenas há que observar 
limites de deformabilidade para assegurar boas 
condições de funçionamento. Assim, é corrente a 
limitação de rotações máximas a valores de 0,5 
a 1 grau sexagesimal em torres e antenas de 
cerca de 100 metros de altura. 


2 — SOLICITAÇÕES DEVIDAS AO VENTO 


No estudo da estabilidade aerodinâmica deste 
tipo de estruturas, há que considerar quatro tipos 
de solicitações : 


— uma solicitação em regime estático, depen- 
dente da velocidade média do vento, da altura e 
da forma da estrutura; 


— uma solicitação em regime dinâmico, cor- 
respondente à turbulência, actuante na direcção 
das rajadas; 


— uma solicitação em regime dinâmico, 
periódica, correspondente à formação de turbilhões 
de Karman, actuante na direcção perpendicular 
à das rajadas, que se manifesta especialmente se 
a estrutura tem a forma circular em planta : 


—e finalmente uma solicitação periódica 
associada ao desenvolvimento dos turbilhões de 
Karman e com frequência dupla destes, que tende 
a ovalizar as secções circulares. 


A quantificação de cada um destes tipos de 
solicitação faz-se do seguinte modo: 


2.1 — Solicitação em regime estático 


Corresponde a uma pressão horizontal actuante 
em qualquer direcção. Segundo o Regulamento 
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de Solicitações em Edifícios e Pontes [1] é da 
forma 
Es 


Re O 6 
16 


res o BB 


* -— 1 
sendo v a velocidade do vento em mseg |, 


c o coeficiente de forma, positivo ou negativo 
conforme se trata de pressão ou sucção, e p a 
pressão em kgfm”*. Nos artigos 49.º a 53.º 
e em anexo o RSEP indica valores para o coefi- 
ciente de forma de diversos tipos de estruturas 
e estabelece valores da velocidade máxima em 
função da exposição das regiões e da altura da 
construção acima do solo. Estes valores resultam 
já da aplicação de critérios estatísticos para 
avaliar a velocidade do vento e correspondem 
a uma representação simplificada dos perfis de 
velocidade de Davenport [2]. 

Segundo o artigo 53.º do RSEP, «para o esta- 
belecimento de coeficientes de forma relativos 
a casos não tratados no anexo, poder-se-ão uti- 
lizar resultados experimentais fidedignos ou 
tabelas geralmente aceites, tais como tabelas 
regulamentares de outros países. 

No caso de construções muito importantes e de 
forma não normal, recomenda-se a utilização 
de ensaios em túnel aerodinâmico». 

Entre as normas regulamentares com maior 
volume de informação sobre coeficientes de forma 
situam-se as normas francesas [3]. Na biblio- 
grafia específica sobre Aerodinâmica encontra-se 
também grande número de resultados de deter- 
minações experimentais do coeficiente de forma |4). 


2.2 — Turbulência longitudinal (na direcção 
das rajadas) 


Os efeitos de turbulência longitudinal, que 
tem sido estudada especialmente mediante a 
realização de ensaios em túnel aerodinâmico são 
quantificados em diversos regulamentos [3,5] sob 
a forma de coeficientes dinâmicos designados por 
coeficientes de rajada, por meio dos quais se 
majora a pressão estática devida à velocidade 
média. 

A quantificação mais perfeita do coeficiente 
de rajada é possivelmente a proposta por 
Davenport |2] e adoptada no regulamento dina- 
marquês de 1967, segundo a qual a pressão 
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devida ao vento, em qualquer ponto de uma 
estrutura, à cota h, é dada por: 


pp causas a 1) 


em que p é a pressão média sobre a estrutura, 
determinada por exemplo a partir de ensaios em 
túnel aerodinâmico tomando em consideração 
a distribuição estatística da velocidade média das 
rajadas, e em que G é o coeficiente de rajada 
dado por 


G=1+grVB+R.... (3) 


sendo : 


g — coeficiente de ponta, função da frequência 
própria fundamental da estrutura, f, e da 
duração T considerada para definir a 
velocidade média (fig. 1); 


| 
dy 
ak] 


0 


Coeficiente de ponta, « 


| f - freq propria fundamental, HZ 


T= duração considerada para 
definir a velocidade media seg 


o 


Fig. 1 — Coeficiente de ponta, g 


r — coeficiente de rugosidade, função da natu- 
reza do terreno e da altura h (fig. 2); 

B — coeficiente de turbulência devida à estru- 
tura, função da altura desta (fig. 3); 

R — coeficiente de ressonância, dado por 
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Fig. 3 — Coeficiente de turbulência da estrutura, B 


em que: 


S — coeficiente função de h/b (relação altura/ 
[dimensão no sentido das rajadas), da 
frequência f e da velocidade média do 
vento v, fig. 4); 

F — coeficiente de energia das rajadas, função 


de F/v, (fig. 5); 
n — amortecimento relativo da estrutura. 
Este critério leva portanto à determinação das 


forças estáticas correspondentes à máxima defor- 
mação. 


2.3 — Turbilhões de Karman 
Sendo R, o número de Reynolds dado por : 


V 
Re=—... 
V 


é a MB) 


TO T 11 = [ELA . . . . — (5 
é st 1 Tu Se E me Sis y — viscosidade cinemática do ar (14,8><10 
o 7h as | NERE pot 9 | 
O = Pa À || | a | 
n os jafega por mi ddd 14 m seg ) 
a ” q Sd ne fas 1) 
z Os TS 42111! . 
2 03 e À E Mt TA e ainda: 
E Da RR o 
a DE 
 Q2 eg q " No é fD 
U NEN bs * V= coesa 8 (6) 
S No: PRN de GR hiso 1 É 
N TE Q Mg t 
| N NON N N 4 
se N MA TN Comi 
N NAN N | — Ago) o 
cai ES q N T E sendo : 
00 N X N bo N | . 4 
si x N A dg] V— velocidade crítica do vento, para a qual a 
a N A q bi] N frequência dos turbilhões é igual à fre- 
nl No K quência própria da estrutura, f (Hz); 
002 Ti 
WIN. ] N 
bg] h” A A S — número de Strouhal: parâmetro adimen- 
001 NV Nha q sional, dependente do número de Reynolds 
aa A N N A e da forma da estrutura, cujos valores se 
PS qa EA, N | indicam na fig. 6, extraída de [6] 
0004 4 E N— 
0003 PAS A formação de turbilhões de Karman dá-se 
a gt ' | | sobretudo para valores de R, > 10“ e a força 
N sinusoidal correspondente, perpendicular à direc- 
N ção do vento, tem o valor: 
eia: 7 O2 03 0405 07 10 20 304050 70 100 
th/v 1 2 
[ *h E(t)=— CV Asen2nft...(7) 
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D— diâmetro do cilindro sobre o qual o vento 
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incide 


seg” à temperatura de 15ºC e pressão 
de 760 mm (Hg). 
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A— área projectada perpendicularmente à 


direcção do vento. 


É de not O 
e notar que —pV = — 
SR e (ar 16 
dezas expressas em kgf, m, seg. 


, para as gran- 


10" 10º 10º 10º 10” 

Re 
Fig. 6 — Valores do número de Srouhal S e do coefi- 
ciente de forma €C em função do número de Reynolds, 
Ra, em cilindros 


2.4 — Dvalização 


Como se referiu, este fenômeno é devido tam- 
bém à formação dos turbilhões de Karman. 
Torna-se necessário evitâá-lo, o que em regra é 
fácil, bastando para isso que a frequência própria 
de ovalização se afaste da frequência dupla da 
de formação de turbilhões. 

A frequência fundamental de ovalização de um 
anel cilíndrico é dada por: 


(a) 


meme 0E es um tb 
sendo: 


— E módulo de elasticidade do material do 
anel; 
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— | momento de inércia polar da secção ; 
— massa por unidade de comprimento; 


R — raio médio. 


3— ANÁLISE ESTÁTICA E DINÂMICA DA 
ESTRUTURA 


Quantificadas as solicitações, a análise aerodi- 
nâmica processa-se ao longo das seguintes 
fases: 


3.1 — Idealização do modelo estrutural 


Adopta-se geralmente como modelo estrutural 
uma viga vertical de secção variável encastrada 
ou apoiada na base e com apoios intermédios 
elásticos correspondentes às espias a vários 
níveis, como se indica na fig. 7. 

O esquema (d) é pouco frequente. Nos casos 
em que se aplica, o apoio tem uma constituição 
tal que funciona como articulado até ao fim de 
montagem e aplicação de pré-tensão nas espias 
e como encastrado posteriormente. 

No caso de chaminés de grande altura, geral- 
mente de betão armado ou pré-esforçado, torna- 
-se necessário considerar a deformabilidade da 
fundação, o que pode ser feito por meio dos 
coeficientes de Barkan [7] e particularmente do 
coeficienle de compressão elástica não uniforme 
c,, visto que a rotação da fundação é o movi- 
mento preponderante. Do ponto de vista de 


idealização estrutural, a deformabilidade das fun- 


dações pode esquematizar-se como se indica na 
fig. 8. 

A rigidez dos tramos fictícios AB deve ser a 
correspondente à deformabilidade da fundação. 
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A deformabilidade por esforço transverso é 
desprezável para este tipo de estruturas. Mesmo 
no caso de chaminês metálicas em que a rigidez 


E ia DS e 
DAR, 
o] 

c 


A A 


Fig. 8 — Idealização estrutural da deformabilidade das 
fundações 


EI é muito grande em relação à área da secção, 
a esbeltez da estrutura leva a que a deformabili- 
dade se deva quase exclusivamente à flexão. 


3.2 — Matrizes de massa, deformabilidade e 
rigidez 


Escolhidos os pontos onde se supõe concen- 
trada a massa da estrutura, o processo mais 
simples de caracterizar a deformabilidade e a 
rigidez é calcular a matriz de deformabilidade 
por aplicação do teorema dos trabalhos virtuais, 
utilizando a representação matricial. 

Considerando a estrutura constituída por ele- 
mentos limitados por nós onde se supõem apli- 


cadas as forças, como se indica na fig. 9, define-se a 


5 

ie (5 

Fa SS 

die 4 

Fa (8) 

Sa 3 Me Ma 
= ad = 

Fx 3) / " I Ae N 

.. to Ea E dd 
(5 De — Bad 

F4 nt 

Ebsagis 1 ) ! | 
1) : 

> O 


Fig. 9 — Consola vertical e tramo elementar 


matriz de deformabilidade de cada tramo |f,| 
tal que 


el=[8 cul o 
ou seja: 
(0)=[6] (M) 


em que f,; é a rotação em k devida à aplicação 
do momento unitário em j. 

Agrupem-se as matrizes [f)] muma matriz 
única |f] tal que 


(10) 


sendo n o número de tramos elementares. Defi- 
ne-se a matriz de transformação de solicitações 
aplicadas em esforços, |[b|, dada por: 


Me | + bi E o iu b” | 

e) IM b) | (E 
Moe | 
Moa ; F, 
oa gel | 

| 
Mae k 

Ma | ba go so va BonÀ 

(11) 

ou seja: 


(Mp= by) (E) 


em que b, é o momento flector em i devido à 
força unitária aplicada em j. 

Demonstra-se facilmente [8 a 10, por aplica- 
ção do teorema dos trabalhos virtuais, que a 
matriz de deformabilidade da estrutura é dada por 


[4] =[b]" LE) Ib] . 0. co (12) 
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Em estruturas isostáticas a determinação de [b] 
não tem qualquer dificuldade. Para estruturas 
hiperestáticas a análise matricial permite calcular 
[b] com relativa facilidade [8 a 10]. No entanto, 
para estruturas do tipo considerado neste capítulo, 
o processô mais simples de determinar [4] quando 
as estruturas sejam hiperestáticas, devido desig- 
nadamente à existência de espias a vários níveis, 
consiste em determinar a matriz de deformabili- 
dade [4º] para uma estrutura base isostática, 
obtida suprimindo uma ou mais ordens de espias. 
Seguidamente, por aplicação do método das for- 
ças, calculam-se as forças de interacção devidas 
às espias suplementares, correspondentes à actu- 
ação de forças unitárias aplicadas nos nós. As 
deformações devidas às forças de interacção apli- 
cadas à estrutura base devem finalmente ser 
adicionadas aos termos correspondentes da matriz 
de deformabilidade daquela para obter a matriz 


de deformabilidade total. 
A fig. 10 exemplifica o que acaba de ser 


exposto. 


obtém-se as restantes colunas da matriz de 
deformabilidade da torre espiada [4]. 
Relativamente aos valores dos termos de [f] e 
[b] neste tipo de estruturas, há que considerar 
o seguinte: 
Em geral é: 
[E]= [m+f oco (15) 


sendo [f]y a deformabilidade por flexão e [f]r 
a deformabilidade devida a esforço transverso. 
Em tramos de secção constante com o vão 1, 


vem 
| Ê | 
GELLI Z 


l 1 —1'| 
dE E 
ir EGA | 1/1 | 
designando por El a rigidez de flexão e por KGA 
a rigidez por esforço transverso, sendo K o coe- 
ficiente de deformação por esforço transverso. 


Para secções anelares de parede fina, K é dado 
por: 


[fly = 
o: « (16) 


4 

A Fa- força de interacção 
devida as espias 

2 

1 


Fig. 10 — Torre espiada e estrutura base 


Considere-se a matriz de deformabilidade da 
estrutura base [A”], A 1.2 coluna de [4º] corres- 
ponde às deformações para uma força unitária 
aplicada em 1. Sob a actuação desta força 
desenvolve-se na espia, em 3, uma reacção F, tal 


que: | 
F, * Aº:s + E, “ k e A ai E (13) 


designando por k a rigidez das espias. Da 
equação 13) calcula-se F.. A deformação da 


od 


torre espiada para F,=1 é então: 
o AO 

S=4 ,—F;. ds 
o A] 

And Es. Ag 
O AO 

Ap =4 ,—F,. A 5 


o o 
A =8 4 —E,. Aa 


- - (14) 


Procedendo análogamente em relação à actua- 
ção de solicitações unitárias em 2, 3 e 4, 
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ala vi) 
4+ 37 


K = (17) 


sendo » o coeficiente de Poisson. 


Em geral pode desprezar-se a deformabilidade 
por esforço transverso. 

Em tramos de secção variável, considerando 
cada tramo constituído por três troços, i, c, j, de 
secção constante, será: 


(0,94, 41,24, — (0,12,40,84 + 


1 = 0,1 x) =| 0,1 2) 
O 
[EJm E o 
SEL] (014 +0,8% + (-0,19+-1,20, + 
+ 0,1 2) E 0,9 x) 
(18) 
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a | 


referindo i, c, j ao troço da esquerda, central e 
da direita e 1, e El, a um dos troços, tomado 
como réferência. 

[ f | continua a ser dado por (16), sendo agora : 


a 4 4,4 ) 
KGA 6 À KGA, KGA,  KGA, 


- o « (19) 


Em qualquer dos casos [f] expressa por (11) 
é susceptível de abaixamento de grau, visto que 
a continuidade exige que a rotação de duas extre- 
midades que convergem num nó seja a mesma. 
Em virtude desta condição simplifica-se [f] supri- 
mindo a última linha e a última coluna (visto 
na extremidade não se desenvolverem em geral 
momentos) e adicionando os termos correspon- 
dentes das linhas 2 e 3) 4 e5,6e 7, etc., pas- 
sando a soma a constituir respectivamente 
as linhas 2, 3, 4, etc. da matriz simplificada. 
Assim por exemplo, a matriz f: 


I Za 1Jà O O 9 0 
1,3 30 0 O O 0 
l Oo 03,923 0 0 
[É] == —— 
6 EI O q BESSA BU 0 
O 0 007,7 4,7 
O DO Gas rs: 


correspondente a uma viga com três tramos 
encastrada na extremidade a que correspondem 
a 1.2 e 2.2 linha e coluna, simplifica-se para: 


| 23 1d O 
SF 1,3 6,9 2,3 
6 

| 23131 4,7 | 


sendo 6,9=3,0 43,9 e 131=54+4-7,7. 


Relativamente à matriz de transformação 
de solicitações em esforços [b], desde que não 
haja momentos localizados nos nós é sempre 
Mie = M, pa OU seja, são iguais os momentos 
das extremidades dos tramos que definem 
um nó, pelo que se podem suprimir as linhas 
de elementos todos nulos e de elementos repeti- 
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dos, baixando o grau de [b] para o número 
de nós. Assim por exemplo, para uma consola 
vertical de três tramos iguais actuada por forças 
horizontais nos nós (fig. 11), em que [b] é defi- 
nida por: 


IM. a 

[Mi 

M, [ 
Bl =] 

IM. | me 

[Mi -| 

(Msa) | 


Fig, 11 


fica simplesmente, em virtude de ser 


FARM ERR 
M, ;-=1I/012 F, 
M,| Joo1l [E] 


No que se refere à rigidez de flexão, expressa 
por [f]!, é de notar que quando a estrutura real 
é triangulada aquela grandeza depende das carac- 
terísticas e da geometria dos elementos que a 
constituem bem como do tipo de ligação entre 
eles. Embora seja corrente admitir uma solidari- 
zação perfeita entre montantes diagonais e 
travessas e calcular o momento de inércia da 
secção composta para essa hipótese, deve-se ter 
em atenção que a rigidez assim obtida é superior 
à real, cuja determinação pode ser feita com 
rigor recorrendo a técnicas analíticas ou experi- 
mentais bem conhecidas e de utilização expedita. 


3.3 — Rigidez das espias 


Há que ter em conta, quando se torna neces- 
sário quantificar a deformabilidade de uma espia 
com as características indicadas na fig. 12, as 
características de comportamento deste tipo de 
elemento estrutural. 

Assim, as equações que relacionam os pará- 
metros referidos na fig. 12 são como se sabe, 
para a hipótese de fios inextensíveis, as seguintes. 


TECNICA Nº 348 


km, ESA 
mi fe A 


* CENTRAL TÉRMICA DA REFINARIA DO NORTE DA SACOR PARA A QUAL À 
MAGUE, EM CONSÓRCIO COM A BROWN BOVERI & CIE. , FORNECEU: 


-4 geradores de vapor de 70t/h cada (produzindo vapor a 68 bar e 455ºC, e 
queimando fuel-oil e gás vindos da refinaria). 
(projecto Mague-Licença Foster Wheeler ) 


=2 turbo-grupos de 14 MW e 12 MW, respectivamente de condensação e cor- 


rapressão. 
(Colaboração Brown Boveri & Cie.) 


= 10 permutadores de calor; 4 torres de destilação e 17 conjuntos de 3 reservo- 
vatórios de alta pressão para hidrogénio. 
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SIDERURGIA NACIONAL 
Laminagem de chapa a frio, 
Tubagens,. 


ALTERNADORES (de potência superior a 2000 kVA), COMPENSADORES 
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L - Desenvolvimento 


A - Área 


E - Modulo de elasticidade 


a, 
H 


(1,60 18x10º Kgfem"2 para cabos de aço) 


p - Peso por unidade de comprimento 


[8 
Cos 


- Carga vertical por unidade de vão 


H - Componente horizontal do esforço axial, 


(constante) 


Fig. 12 — Características das espias inclinadas 


a) ue [AD qe É cold 

| ses ara 
L=1(1+ apo 

? | 21 


8 fé cos" 2 ) 
À 5º 


É usual adoptar-se para tensão inicial na espia o 
valor correspondente à flecha máxima des/40 a s'80. 

Quando a extremidade superior da espia 
é sujeita a um deslocamento horizontal à, a varia- 
ção das forças H e T é tal que a deformabilidade 
resultante À, é dada por [11]: 
3 á p? 6 
da T* 


À = ae AB 


ton so 
i AH EA cos" « 


A primeira parcela de À, resulta do esforço axial 
devido ao alongamento da espia considerada como 
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E 


barra recta sem peso (fig. 13) e é sensivelmente 
constante, visto que a variação de «a é desprezável. 


Fig. 13 — Alongamento da espia considerada como barra 
recta sem peso 


Da fig. 13 resulta que: 


RD AO E ge E 
AH AT cos As EA cos“ 
| | 
a er é aa) 
EA cosa 


A segunda parcela de (21) traduz a não-linea- 
ridade do comportamento da espia e varia rapi- 
damente com a tensão T. Para garantir que 
a espia oponha uma determinada resistência 
a deslocamentos transversais impostos à extre- 
midade superior há pois que estabelecer criterio- 
samente a pré-tensão inicial T, a instalar. 

Como a rigidez da espia aumenta rápidamente 
com o deslocamento ), na análise estática ou 
dinâmica de um mastro espiado a vários níveis 
há que estabelecer critérios que permitam tratar 
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a espia como um apoio elástico, dentro de uma 
aproximação razoável. A fig. 14 ilustra a aplica- 
ção deste tipo de critérios. 

No caso de espiamento a um único nível os 
esforços no mastro são independentes da defor- 
mabilidade das espias, pelo que as considerações 


da espia H 


ad 


Reacção 


(fig. 15) cujo estudo específico se torna frequen- 
temente necessário. 

Para quantificar a rigidez de torção são válidas 
as considerações feitas no decorrer deste pará- 
grafo relativas à rigidez de cada espia. Há que 
compor a rigidez de todas as espias a um dado 


comportamento rea! (eq 21) 


tg $ = rigidez da espia linearizada 


deslocamento 5 


deslocamento 
maximo previsto. 


Fig. 14 — Critério de linearização da espia 


feitas só apresentam interesse prático desde que 
seja necessário limitar as deformações do mastro, 
o que é corrente, como se referiu em 1. 

Embora seja possível em teoria realizar por 
intermédio de espias apoios com uma rigidez 
préviamente estabelecida, na prática há que ter 
em conta que o acréscimo de rigidez é conseguido : 
ou à custa do aumento da secção dos cabos, o 
que é condicionado pela existência no mercado ; 
ou diminuindo a inclinação x, o que é condicio- 
nado pelo espaço disponível; ou aumentando a 
pré-tensão, o que implica espias mais fortes, 
maiores esforços nas respectivas fundações e com- 
pressões instaladas no mastro, 

A rigidez do espiamento a um determinado 
nível depende do número e da disposição das 
respectivas espias; existem vários métodos para 
determinar a rigidez de conjunto em relação a 
deslocamentos numa determinada direcção (!) 
[12]. 

É de notar ainda que na ligação das espias ao 
mastro se utilizam correntemente dispositivos 
destinados a garantir a necessária rigidez de torção 


(1) Na publicação referida [12), analisa-se também o 
estudo dos efeitos térmicos no mastro e na espia, que não 
se aborda aqui por transcender o âmbito deste artigo. 
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nível para obter a rigidez de torção conjunta a 
esse nível [12]. 


mastro 


estrutura de ligação 


ALÇADO 


PLANTA 


Fig. 15 — Disposição característica de seis espias para 
garantir rigidez de torção 
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Quando o deslizamento transversal imposto à 
extremidade da espia é da forma 


3 sen wt 
sendo » geralmente a frequência própria de estru- 


tura, a rigidez da espia em regime dinâmico, K, 
é dada por [63]: 


1 N OE 
k=,-=K,| 1-4 =t | . (23) 
d . Goó(U)=1 
em que: 
1 EA cos a 
q 22 
| A, a 
r — sena E 
2 ps EK, 
MR 
9 j 
. = 
G =—s— 
s pk, 


CR a TT 
tg 
2 


EO = D* | ===> 


= 


2 


sendo g a aceleração da gravidade e tendo sido 
os restantes simbolos definidos anteriormente. 

Na prática, a redução da rigidez devida ao 
regime dinâmico expressa pelas equações (23) e 
(24) só tem importância em geral para espias com 
mais de 100 metros de comprimento. 


3.4 — Deslocamentos e esforços 


Efectuada a análise estrutural para os diferentes 
tipos de solicitações referidos em 2., verifica-se 
que as características globais de comportamento 
no que diz respeito a deslocamentos e esforços 
são as seguintes: 

Em chaminés metálicas cilíndricas, tronco-=cóni- 
cas ou mistas (tronco-cônicas na base e cilin- 


TECNICA Nº 395 


dricas na restante altura), a frequência própria 
fundamental é geralmente inferior a 1 Hz e a 
velocidade critica do vento para a formação de 
turbilhões de Karman é da ordem de 20 a 30 km h. 
Como o amortecimento Y é muito pequeno, 
o coeficiente de amplificação em ressonância 
M==1 27 atinge facilmente valores superiores 
a 100, o que leva a deslocamentos e esforços 
elevadíssimos, que excedem a capacidade resis- 
tente da estrutura. O efeito Karman é então de 
primordial importância e torna-se indispensável 
evitá-lo, utilizando os processos a que se fará 
referência em 4, 

A actuação do vento excepcional neste tipo de 
estruturas é feita principalmente em regime está- 
tico e os esforços correspondentes são em geral 
facilmente absorvidos por espessuras de chapa 
pequenas. 

É de salientar que a formação dos turbilhões 
de Karman nunca é simultânea com a ocorrência 
de vento excepcional, dados os baixos valores 
da velocidade crítica. 

Em chaminés de betão armado o valor do 
amortecimento é consideravelmente superior ao 
das chaminés metálicas, pelo que o efeito de 
Karman, embora frequente, pode não levar 
a deslocamentos e esforços exagerados. Como 
critério de julgamento sobre se deve ou não ser 
evitado pode empregar-se o que consiste em 
comparar os seus efeitos com os do vento habi- 
tual, desde que não sejam simultâneos. Se os não 
excederem, pode permitir-se a sua formação. 

O dimensionamento de chaminés de betão 
armado e pré-esforçado põe um certo número 
de problemas específicos, tratados por exemplo 
em [13]. 

Em torres e antenas é predominante a actuação 
do vento excepcional, sobretudo pela necessidade 
de limitar as deformações. Os efeitos de turbu- 
lência são importantes, conforme se conclui da 
aplicação de métodos analíticos ou experimentais 
ou de critérios simplificados para a quantificação 
do coeficiente de rajada, referidos em 2.2. 

No que se refere ao dimensionamento das 
espias, este é condicionado pela rigidez que se 
pretende conferir por seu intermédio à estrutura, 
pela respectiva disposição e número, pela pré- 
“tensão inicial que é normalmente estabelecida de 
modo que a flecha máxima (eg. 20) tenha valores 
compreendido entre 1/40 e 1/80 do comprimento 
da espia e ainda eventualmente pela necessidade 
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de se considerarem as características dinâmicas 
das espias. Estas materializam cordas vibrantes 
cuja frequência própria e dimensões podem ser 
tais que levem à formação de turbilhões de Kar- 
man. Tem-se verificado este fenómeno sobretudo 
em zonas onde a formação de gelo e de neve 
sobre as espias dá origem a acréscimos conside- 
ráveis na massa e no diâmetro aparente destas. 


4 — MÉTODOS PARA ATENUAR OU ELIMI- 
NAR A FORMAÇÃO DE TURBILHÕES 
DE KARMAN 


Consegue-se evitar os efeitos dos turbilhões de 
Karman ou por alteração das características dinã- 
micas da estrutura ou por meio de dispositivos 
que impeçam a formação de turbilhões 

A alteração das características dinâmicas con- 
sistirá em geral no reforço das secções ou do 
espiamento de modo a aumentar considerável- 
mente a frequência própria da estrutura, levando 
a que a velocidade crítica passe a corresponder 
a uma probilidade de ocorrência diminuta. Em 
chaminés metálicas este processo não é normal- 
mente viável, pois tornaria necessário aumentar 
a frequência própria para mais de 3 vezes do seu 
valor inicial, o que significa aumentar a rigidez 
da estrutura mais de 9 vezes. 

O processo mais prático para evitar ou atenuar 
o formação de turbilhões consiste em dispor 
sobre a estrutura elementos isolados ou ligados 
formando hélices [14 a 16]. O estudo da eficiência 
de vários tipos de hélice tem sido feito, com carác- 
ter sistemático, no National Physical Laboratory 
em Inglaterra, tendo sido publicadas, no número 
de Outubro de 1966 da revista «Engineering 
News Record», as seguintes especificações a adop- 
tar em chaminés: 


passo dos hélices: 5 D (D — diâmetro exterior 
da chaminé); 

altura dos hélices: terço superior da chaminé; 

número de hélices: 3, igualmente espaçadas 
em planta; 

altura do perfil dos hélices: 0,09 D. 


A fig. 16 exemplifica estas especificações 


Julga-se que este dispositivo foi aplicado pela 
primeira vez em Portugal em 1966 numa chaminé 
de um complexo industrial próximo da Figueira 
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5 D 


PLANTA 


dg 


Fig. 16 — Características de hélices para dissipação dos 
turbilhões de Karman 


da Foz (fig. 17). De então para cá tem sido utili- 
zado em diversos locais. 


Fig. 17 — Chaminé protegida por hélices 
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NOTICIÁRIO 


2.º Congresso da Sociedade Internacional 
de Mecânica das Rochas 


e cem 
- — 


Este Congresso decorrerá em Belgrado, de 21 
a 26 de Setembro de 1970. Na sua agenda 
incluem-se os seguintes temas : 


1 — Propriedades intrínsecas dos maciços rocho- 
sos. 

2 — Deformabilidade dos maciços rochosos. 

3 — Existências mecânicas dos maciços rochosos. 

4 — Trabalhos subterrâneos. 

5 — Cominuição. 

6 — Estabilização dos maciços rochosos. 

7 — Estabilidade dos taludes naturais e de esca- 
vações permanentes e temporárias. 

8 — Comportamento dos maciços rochosos 
quando sujeitos a determinados tipos de 
solicitação. 


Durante este Congresso estará patente uma 
exposição cientifica e comercial de material de 
campo e laboratório que se utiliza no domínio 
da mecânica das roçhas. 

Toda a correspondência referente a este Con- 
gresso deverá ser dirigida a: 


Secretarijat 2. Kongresa Medjunarodnog 
Drustva Za Mehaniku Stena — Institut Za 
Vodoprivredu «Jaroslav Cerni», Bulevar Voj- 


vode Misica 43, Beograd, Jugoslavija 


Conferência sobre materiais dielectricos, 
medidas e aplicações 


e— = 


A realizar de 20 a 24 de Julho de 1970, na 
Universidade de Lancaster, Baibrigg, Lancaster, 
esta conferência realiza-se sob os auspícios de 


Institution of Electrical Engineers, 

Institute of Physics, 

Physical Society, 

Institute of Electrical and Electronics Engineers. 


e interessa a quem se dedique à pesquiza e apli- 


cações de dieléctricos sólidos e líquidos. 
Para outros esclarecimentos, dirigir correspo- 
dência a: 


Mrs. Ann Cook 
Manager, Conference Department, 
The Institution of Electrical Engineers, 


Savoy Place, London WC2R OBL 
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«O Mundo à água e nós» 


exposição internacional na poluição de água e tratamento 
requerido, de 2 a 9 de Setembro, de 1970 em Johkopihg, 
Suécia 

Os problemas da conservação da água e da 
natureza são hoje discutidos em todo o mundo 
e a protecção das águas é frequentemente o tema 
central nestes campos. 

As actividades intensivas de informação e a 
crescente cooperação entre países para resolver 
problemas técnicos e práticos é hoje imperativa. 
O projecto «O Mundo, a água e nós» é um 
importante passo nesta direcção. O organizador 
é Elmia Ab Jôhkóping na Suécia, em cooperação 
com o Secretariado de Conservação da Natureza 
Sueca, em Estocolmo e o Comité Nacional Sueco 
para o Ano Europeu de Conservação das Aguas. 

«O Mundo, a águas e nós» é uma figura inter- 
nacional especializada para a conservação das 
águas com uma conferência internacional e expo- 
sições de ideias, onde as organizações e as com- 
panhias apresentarão os seus produtos e os seus 
acabamentos dentro do campo da conservação 
de águas. 

«O Mundo a água e nós» inclue; 

Exposição de ideias tocando os seguintes tópi- 
cos : 


1.º — A água e um problema para a raça 
humana: 

2.º — A utilização da água 

3.º — A comunidade e a conservação da água 

4.º — Os lagos de Watterh Wãherh — Mailareh, 
descrição e comparação entre os polui- 
dores dos três lagos 

5.º — A necessidade para experiências, treino 
e informação geral. 


O projecto e a exposição de produtos: 


1.º — Purificação biológica e química 

2.º — Tratamentos através da Engenharia de 
(plantas de purificação, superfície da água 
e distribuição) 

3.º — Processos e contrôle químico 

4.º — Irrigação, medidas laboratoriais e técnica 
de restauração 

5.º — Petróleo e sua pesquisa. 


Conferências 


As conferências serão de dois tipos onde se 
tratará da preservação da água com orientação 
no mar Báltico, dos estudos nos grandes lagos 
bem como a discussão dos problemas expostos 
anteriormente. 

Mais esclarecimentos: World Water And We 
Sweden — Box 6066 5-55006 Jônkôping 6 
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A indústria de 
embalagem 
aprecia-nos 
porque poupamos... 
poupamos 

o seu dinheiro! 


A MOBIL pode demonstrar-lhe como 
gastar menos produto no revestimento de emba- 
lagens. 
Desde o cartão canelado às embalagens mais sofisticadas, 
os MOBIL-KOTE podem proporcionar as características de 
protecção e aparência exigidas em cada caso. 

Os Engenheiros da MOBIL mostrar-lhe-ão como econo- 
mizar dinheiro na sua Indústria, aplicando a menor quantidade 
possível de MOBIL-KOTE. 


Mobil 


vendemos mais vendendo menos 


AMobil = 
mostra-lhe 


inumeras maneiras 


es —— — 
ss — qi 


Por ano poupamos à 
industria vários milhares 
de escudos. 

Uma grande parte 

deste dinheiro 

sai das nossas algibeiras. 


vendemos ma:s vendendo menos 


Os Mobiltac 
ventam-se 
tanto em serviço, 
que dificilmente 
poderemos 


vendé-los muito. = 


vendemos mois vendendo menos 


Mobil 


vendemos mais vendendo menos 


e 


Este simples ensaio 
mostra como reduzir 

os custos 

de manutenção 

dos sistemas hidráulicos. 


o mo O a e am 
na O e — 


vendemos mais «sendendo menos 


A industria de 
embalagem 
aprecia-nos 

porque poupamos... 
poupamos 

o seu dinheiro! 


vendemos mas vendendo menos 


Se é raro apresentarmos 
o preço mais baixo, 
porque razão vendemos 
tanto óleo para 

grandes turbinas? 


vendemos ma:s vendendo menos 


C. D. U. 627.825:627.824 


O MÉTODO DE AJUSTAMENTO COMPLETO 
PARA O CÁLCULO DE BARRAGENS ABÓBADA 


(Conclusão) 


7 — CÁLCULO DE ARCOS 


Os princípios básicos usados para o cálculo 
matricial dos arcos, de modo a determinar as 
forças internas e os deslocamentos devidos a 
cargas unitárias nas aduelas, já foram explicadas 
noutro lugar**. A seguir dá-se um sumário da 
teoria que ultimamente foi desenvolvida em resul- 
tado dessas investigações já referidas. Na Fig. 5 
indicam-se os dados geométricos usados no cál- 
culo de um semi-arco, através de um programa de 
computador electrónico. 


Fig. 5 — Notações para um semi-arco 


A definição de cada aduela de um arco é feita 
através dos seguintes elementos; e, — espessura 
na secção i; R,;—raio da face de montante da 
aduela; ó,, — ângulo entre o eixo da aduela ou 
secção i e o fecho f do arco; e x, e y; as coórde- 
nadas da secção i em relação ao fecho. O arco é 
dividido em i aduelas numeradas de O (zero) 
até f. 
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As características da fundação, necessárias para 
o cálculo, são expressas pela matriz [C,]| (ver 
Ref. 6 e 15). 


(38) 


A partir dos dados geométricos define-se a 
matriz [B, 1] que transforma os esforços na secção 
j, à direita da secção i, em esforços na secção i, 
EÉ 


[B,] = | 0 cos 6; sen 6, (39) 


O — sen ur cos ur 


sendo M os momentos de flexão — positivos no 
sentido dos ponteiros do relógio; N os impulsos 
(tangenciais em relação à linha do centro de 
gravidade da secção) — positivos quando dirigidos 
no sentido do encontro; e Q as forças de corte 
(tangenciais à secção ou radiais ao arco) — que 
são positivas quando tenham o sentido montante- 
-jusante. 

A matriz transposta [B,]* transforma os 
deslocamentos nas secções i em deslocamentos 
nas secções |. 
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Para cada aduela escrevem-se também as 
seguintes matrizes (de três linhas e colunas): 


Matriz [C] — que exprime os deslocamentos da 
secção limite da aduela, devidos 
a cargas unitárias concentradas, 
admitindo a aduela livre nesta 
secção e fixa na secção oposta; 

Matriz [D] — que exprime os deslocamentos nas 
condições anteriores mas devidos 
a cargas unitárias distribuídas 
actuando na aduela; 

Matriz [E] — que exprime os esforços na secção 
à esquerda da aduela, devidas às 
cargas unitárias distribuídas na 
aduela, considerada livre à direita 
e fixa à esquerda. 


Em publicações anteriores descreveram-se mais 
completamente as expressões de cada um dos 
elementos das matrizes anteriores. Para cada 
secção i pode ser deduzida a seguinte matriz: 


(ad = 5 Bd [Co (Bdt| (87 


(40) 


que exprime o deslocamento do fecho do arco 
devido à carga concentrada unitária na secção i. 
Os deslocamentos [h,] do fecho do arco 


devidos à carga uniformemente distribuída na 
aduela (i—1l,i), são: 


[hal=1 A, ]><LE1+4[B;I ID, (41) 


UMA UNIDADE | 
DE CARGA 


LINHA DE CENTRO 


DC ARCO 


Fig. 6 — Distribuição das cargas unitárias num arco 


A fig. 6 indica uma parte do arco represen- 
tando a posição das consolas e a forma das 
cargas triangulares unitárias tendo os seus vérti- 
ces nos pontos o, m, j, n e f que são pontos de 
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intersecção do arco com as consolas. Para cada 
aduela i podem-se obter dois coeficientes 1, e 1; 
que dependem do comprimento da aduela e das 
posições das consolas, à direita e à esquerda 
dessa aduela. Tendo em consideração a Eqg. 41, 
os deslocamentos [h,] no fecho do arco devido 
às cargas triangulares unitárias (U,| podem ser 
determinados pela expressão seguinte : 


j n 
[hJ=Sixlhl+ S xlha] (42) 
i=m I=) +] 

Uma vez calculados os deslocamentos [h,] no 
fecho do arco, para cada carga triangular unitária, 
é com facilidade que se determinam as incógnitas 
no fecho (M, N e T) para cada uma das cargas 
consideradas. Usando a matriz de transformação 
[BJ], as incógnitas em cada secção i para cada 
carga unitária (U |, são obtidas pela expressão : 


M M 
N = [B,] x N + 2 LxI[B;i-n]> 
V Vl t=m 


x[El+ 5 Lx [B,q-]> [E] 
i=j+1 
| (43) 


Para cada tipo de carga unitária considerada 
(radial, tangencial, de torsão ou de flexão) o con- 
junto de forças desconhecidas é determinado 
para cada secção j. Uma vez determinadas essas 
forças, os coeficientes k , que entram no sistema 
de Eg. 36, podem ser determinados pela expres- 
são seguinte: 


o MI 
[k,4] = [B,;] ne ain e 
V Ju 


(44) 


2 h5 bd 


Ii j-l 


+ 21x [h]+ 


I=m 


8 — CÁLCULO DAS CONSOLAS 


O cálculo das consolas é muito parecido com 
o cálculo dos arcos, mas considera as consolas 
como elemento torcidos pelas razões atrás expos- 
tas. Cada consola (Fig. 7) é normalmente consi- 
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derada dividida em secções j correspondendo 
as suas intersecções com os arcos e às secções 
intermédias. A definição dessas secções é comple- 
tamente obtida pelos dados seguintes: 


Fig. 7 — Notações para uma consola 


da face de montante à superfície de referência ; 


y,—— distância vertical do coroamento até 
à secção; R uj — Taio horizontal da face de mon- 
tante na secção j; é, — ângulo da direcção radial 
da consola com o eixo de simetria (ou com 
a consola central) da barragem ao nível da sec- 
ção j. 


e, — espessura na secção j; m, = distância 


Para iniciar o cálculo determinam-se não 
só as forças internas M, e N, (momento flector 
e impulso), devidas à subpressão na base, mas 
também os coeficientes de deformação da fun- 
dação (2,,%+ 35, Bo)”. 

As consolas são normalmente calculadas para 
o seu peso próprio e para a componente verti- 
cal da pressão hidrostática, associadas com a sub- 
pressão na base. No que se refere à componente 
horizontal da pressão hidrostática, assim como 
às forças internas de interacção com os arcos 
(forças verticais e flexões) são consideradas car- 
gas triangulares unitárias. Elas devem ter um valor 
unitário em cada ponto de cruzamento com o arco 
e devem ser iguais a zero nos pontos de cruza- 
mentos com os arcos acima e abaixo do anterior. 

Determinam-se, em cada secção j, para cada 
tipo de carga considerada, as seguintes forças 
internas : 


N, — impulso (vertical) 


Q, — esforço transverso radial (horizontal) e 
M, — momento flector (de eixo horizontal) 
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Para cada secção, tomando em consideração 
os dados anteriormente descritos, calculam-se, 
também, as seguintes características : 


S; — a área da secção; 
1, — momento de inércia da secção ; 
d;, — distância do centro de gravidade à face 


de montante 


Para calcular os deslocamentos k,, devidos 
a cada carga, usa-se a expressão geral que é 
deduzida fácilmente da teoria da flexão: 


j 

1 pá M , 
+ DO MY a (RE) cos (6, ) + 

n=j--+-1 n 


M, (Y,— Yp-4) 3F, + Fo 


— E,) + 
I, 2 


+M, 8, (E — E) +LIç(Y,—Y) cos (6, — 6) + 


1H b-1 H 
+= -Y% pit > gar Canal 
- 9; n=j+1 Un 
H, 
E Lg ei! (45) 


5» 


Nesta expressão b designa a secção da base 
da consola e as funções F, He L são respecti- 
vamente : 


para deslocamentos verticais z—d, E = ns 


para deslocamentos radiais y |, 


K 05 

Eee = —— ' Q 

G e E ds 
considerando v=0,22eK ==1,5 


Também se determinam, em cada ponto de 
cruzamento com os arcos, rotações de eixo ver- 
tical devidas às cargas de torção triangulares do 
mesmo tipo que as indicadas anteriormente. 
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9 — COMPUTADORES USADOS 


O cálculo anteriormente descrito foi inicial- 
mente programado para um computador Elliott 
803-B. 

Em geral são usados cinco programas suces- 
sivos : 


1) Características geométricas dos arcos e das 
consolas obtidas através da definição matemática 
das formas? e divisão destes elementos em aduelas; 

2) Cálculo de arcos para três tipos de cargas 
unitárias (radiais, tangenciais e de flexão de eixo 
vertical); 

3) Cálculo das consolas para quatro tipos de 
cargas unitárias (radiais, verticais, flexões hori- 
zontais e torsões de eixo vertical); 

4) Estabelecimento e resolução do sistema de 
equações ; 

5) Cálculo dos deslocamentos, tensões normais 
e tangenciais nos arcos e nas consolas e cálculo 
das tensões principais nos paramentos. 

Estes programas estão agora a ser usados no 
novo computador NCR-4100, depois de terem 
sido traduzidos em ALGOL. Este computador 
poderá realizar o cálculo completo numa única 
operação. 


10 — COMPARAÇÃO DOS RESULTADOS 
OBTIDOS POR AJUSTAMENTO COM- 
PLETO COM OS RESULTADOS OBTI- 
DOS POR OUTROS MÉTODOS 


Uma das melhores verificações sobre a precisão 
do método de ajustamento completo é a que se 
obtém pela comparação dos resultados obtidos 
com vários métodos apresentados ao Simpósio de 
1968 de Barragens Abóbada”. 

Esta reunião foi realizada na Institution of 
Civil Engeneers, de Londres, para apresentar e 
discutir os resultados de um programa de inves- 
tigação, neste campo, com 10 anos, realizado em 
várias universidades britânicas. Este esforço teve 
o apoio financeiro da CIRIA, «Construction 
Industry Research and Information Association». 

No Simpósio foi apresentada uma comparação 
crítica de resultados dos vários métodos, com 
ensaios de modelos para quatro barragens típicas, 
através de diagramas de resultados. Esses diagra- 
mas de cálculo são agora aqui comparados com 
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os diagramas obtidos pelo método de ajustamento 
completo, para as barragens tipo 1 e tipo 5. 
A barragem tipo 1 é uma barragem cilíndrica, de 
espessura constante, inserida num vale rígido 
trapezoidal. A barragem tipo 5, é uma barragem 
de dupla curvatura, com arcos de espessura cons- 
tante, muito semelhante às formas da barragem 
do Cabril, especialmente à do paramento de 
montante, e também inserida num vale trapezoidal. 
A barragem tipo 1 foi considerada com 30 metros 
de altura e a barragem tipo 5 com 120 metros de 
altura (6 unidades de comprimento) (ver Ref. 5). 
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Fig. 8 — Barragem CIRIA Tipo-r — Comparação dos 
deslocamentos 


As figs. 8 e 9 apresentam os valores dos deslo- 
camentos radiais, tangenciais horizontais e verti- 
cais ao longo do coroamento e do eixo da barra- 
gem (consola central), no caso da barragem 
tipo 1, obtidos pelo método de ajustamento com- 
pleto e por métodos britânicos e, também, os 
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obtidos pelo método «trial load», só com ajusta- 
mento radial, calculados no Bureau of Reclama- 
tion. Como se pode verificar nestas figuras, os 
resultados do método de ajustamento completo 
concordam, muito exactamente, com os obtidos 
com os outros métodos e, em especial, com o 
método de relaxação dinâmica. 

Também se pode concluir que as diferenças 
entre os vários métodos não são significativas, 
com a excepção dos resultados obtidos pelo «trial 
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Fig. 9 — Barragem CIRIA tipo 1 — Comparação dos 
deslocamentos verticais € tangenciais horizontais 


load» radial. Isto indica que se pode considerar 
já realizada a solução matemática do problema 
estrutural das barragens em abóbada e em casca. 

A fig. 10 apresenta as tensões normais, verti- 
cais e horizontais, no eixo da barragem (consola 
central). Também aqui o acordo entre os vários 
métodos é excelente, notando-se uma concordância 
quase perfeita entre o método de ajustamento 
completo e o método da relaxação dinâmica. 

Os únicos valores completos das tensões nos 
paramentos da barragem, que foram apresentados 
durante o Simpósio, foram os obtidos pela relaxa- 
ção dinâmica *”, tendo esse método indicado as 
tensões principais em termos das linhas de igual 
valor da tensão principal máxima e da tensão 
principal mínima em cada ponto dos paramentos. 
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Fig. 10 — Barragem CIRIA tipo x — Comparação de 


tensões horizontais e tensões verticais normais 
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Fig. 11 — Barragem CIRIA tipo 1 — Comparação 
das tensões principais 


A Fig. 11 apresenta a comparação dessas 
linhas de iguais valores das tensões principais 
máximas, obtidas por esse método, em compa- 
ração com os valores das tensões em pontos 
de cruzamento arco-consola obtidas pelo método 
de ajustamento completo. Como se pode verifi- 
car, os resultados têm o acordo perfeito em toda 
a região central da barragem. De facto, nessa 
região, as diferenças entre as tensões principais, 
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Fig. 12 — Barragem CIRIA tipo 5 — Comparação 
dos deslocamentos radiais 


em ambos os métodos, não são superiores, em 
nenhum caso, a uns poucos kg'cm*, Todavia, 
próximo das fundações, há discrepâncias que 
podem ser atribuídas ao carácter de cada método 
e também à dimensão das malhas usadas num e 
noutro método. Há que ter em conta que, junto 
das fundações, o método de ajustamento com- 
pleto admite uma distribuição linear de tensões, 
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mesmo na própria superfície de fundação, 
enquanto que o método de relaxação dinâmica, 
como o método da teoria da elasticidade, chega 
a tensões que teriam valores infinitos nos ângu- 
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Fig. 13 — Carragem CIRIA tipo 5 — Comparação dos 
deslocamentos verticais e tangenciais horizontais 


los entre os paramentos e as fundações, se a 
malha fosse infinitesimal. Para malhas tendo 
dimensões semelhantes às distâncias entre arcos 
e consolas, as tensões deverão ser semelhantes, 
o que, de facto, se verifica. 

No que se refere à barragem do tipo 5, nas 
Figs. 12 e 13 comparam-se os deslocamentos 
obtidos pelos vários métodos apresentados no 
Simpósio *, com os obtidos pelo método de ajus- 
tamento completo. Também aqui o acordo dos 
resultados com a relaxação dinâmica é muito 
melhor do que com qualquer outro método. Pode- 
-se observar que as tensões normais horizontais, 
ao longo do eixo da barragem (consola central), 
e as tensões normais verticais na mesma secção, 
têm um acordo bastante perfeito (Fig. 14). 

As únicas tensões nos paramentos da barragem 
tipo 5, que foram apresentados durante o Sim- 
pósio, foram as obtidas pela relaxação dinâmica, 
e essas tensões foram apresentadas em termos de 
tensões mormais verticais e tensões normais 
horizontais” (Figs. 15 e 16). De igual modo, 
essas tensões são comparadas nessas figuras com 
os valores obtidos pelo método de ajustamento 
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Fig. 14 — Barragem CIRIA tipo 5 — Comparação 
das tensões verticais e horizontais normais 
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Fig 15 — Barragem CIRIA tipo 5 — Comparação 
das tensões normais verticais 
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Fig. 16 — Barragem CIRIA tipo 5 
Comparação das tensões normais horizontais 


completo. Uma vez mais o acordo entre os resul- 
tados é excelente para toda a parte central da 
barragem, ainda que os dois métodos difiram 
bastante nos seus princípios. Finalmente, a 
Fig. 17 apresenta as tensões principais nos para- 
mentos obtidas pelo método de ajustamento 
completo, não tendo sido possível fazer compa- 
ração dessas tensões principais com valores obti- 
dos por qualquer outro método porque, até hoje, 
ainda o não conhecemos que tais valores de ten- 
sões principais tenham sido publicados. 


11 — INFLUÊNCIA DAS DISTORÇÕES DEVI- 
DAS AO CORTE RADIAL E INFLUÊN- 
CIAS DA CURVATURA DA ADUELA 


No Cap. 3 faz-se referência ao facto de que as 
distorções angulares à, e à, poderiam ser despre- 
zadas nos cálculos, especialmente no caso de 
barragens finas. Os cálculos cujos resultados 
foram apresentados no Cap. anterior, para a 
barragem do tipo 5, permitem uma estimativa 
da importância daquela hipótese. Para esse efeito, 
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Fig. 17 — Barragem CIRIA tipo 5 — Tensões principais 
obtidas nos paramentos pelo «Metodo de ajustamento 
completo» 
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Fig. 18 — Barragem CIRIA tipo 5 
Comparação das rotações e distorções angulares 
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calcularam-se os valores das rotações nas secções 


horizontais e verticais, “bia e de Li e também 


d x dy 

os valores das distorções angulares 2, e d,. Estas 
distorções foram obtidas a partir dos valores das 
tensões tangenciais. A Fig. 18 indica os diagra- 
mas daquelas rotações e o diagrama das suas 
somas com as correspondentes distorções angu- 
lares. Como os dois diagramas diferem apenas 
muito ligeiramente, pode concluir-se que o efeito 
das distorções angulares é, em geral, desprezável. 

Também se considerou com interesse determi- 
nar a ordem de grandeza do erro cometido 
quando N,, e N,, são considerados iguais entre 
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Fig. 19 — Barragem CIRIA tipo 5 — Influência da 
curvatura nos momentos de torção e esforços trans- 
versos que actuam nas faces horizontal e vertical das 
aduelas de cruzamento arco-consula tangenciais nos 

paramentos de montante e jusante 
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si (e iguais a T), i.e. quando a curvatura da 
barragem pode ser desprezada nas condições de 
equilíbrio das aduelas (ver Eg. 31 e 9f),e ainda 
quando se consideram M,, eM,, iguais entre 
si (e iguais a M,). Com este propósito foram 
calculados os valores das relações de N ,, 'Nyy € 
de M,, M,, que são apresentados, na Fig. 19, 
Os valores de Mx e de N ,, foram obtidos a partir 
das Eg. 24 e 25. Uma vez que os valores das rela- 
ções atrás referidas são muito próximos da uni- 
dade, pode concluir-se que a hipótese anterior é 
aceitável. Na Fig. 19 podem observar-se três dis- 
crepâncias que necessitam de uma explicação. 
No cruzamento C3-A2 o valor de M,, é muito 
pegueno, o que faz com que a relação com M,, 
seja muito maior do que a unidade. No que res: 
peita às relações entre N,, e Nyy NOS cruza- 
mentos AS—C3 e A6—C2, o erro resulta do 
facto de que as derivadas dos momentos, pró- 
ximo dos extremos dos arcos, são determinadas 
com uma precisão que depende da distância 
entre as consolas. Já se verificou que se se tomar 
uma malha de arcos e consolas mais fina do que 
a usada neste caso, este erro resulta muito mais 
pequeno. 
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Por outro lado ainda, com a finalidade de 
estimar a influência da hipótese de que N,, 
= Nox = T é Mg = Mox = M,, nos resultados 
finais, calcularam-se os valores das tensões de 
corte nas faces horizontais e verticais da aduela. 
Na Fig. 19 mostram-se os valores destas tensões 
tangenciais nos paramentos de montante e de 
jusante da barragem. Pode verificar-se que os 
valores das tensões tangenciais nas faces hori- 
zontais e verticais são muito próximos um do 


outro. Conclui-se daqui que o erro das tensões 
principais, que resulta da hipótese que vimos 
fazendo, é desprezável. Uma vez que a barragem 
tipo 5 tem uma espessura e uma curvatura de 
valor médio, parece que se pode tirar a conclusão 
de que aquela hipótese ainda é aceitável em 


barragens relativamente espessas. 


12 — COMPARAÇÃO COM RESULTADOS 
DOS MODELOS 


Recentemente foram projectadas um número 
relativamente elevado de barragens de arcos de 
três centros?, algumas das quais foram objecto 
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Fig. 20 — Barragem de ATAZAR - Tensões principais obtidas pelo método de ajustamento 
completo e por ensaio em modelo 
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de ensaios estruturais de modelos e foram também 
calculadas pelo método de ajustamento completo. 
Os resultados obtidos para a barragem de Atazar 
(de 135 m de altura) e para a barragem do Funil 
(de 92 m de altura), respectivamente na Espanha 
e no Brasil, foram comparados, noutra publicação, 
e aqui são indicados nas Figs. 20 e 21. Como se 
verifica, os valores máximos das tensões são 
aproximadamente, da mesma ordem de grandeza 
tanto pelo modelo como pelo cálculo, e tanto no 
fecho da barragem como junto das nascenças. 
No caso da barragem do Funil, pode dizer-se 
que o acordo é excelente e que os valores 
máximos das tensões, tanto no fecho dos arcos, 
como ao longo dos limites da casca, são muito 
aproximadamente iguais. 

Comparações semelhantes foram feitas, também 
recentemente, para a barragem de Montanejos 
(de 180 m de altura). Nesta barragem a fundação 
tem um módulo de elasticidade igual ao do 
betão, tendo os modelos estruturais sido realizados 
e ensaiados no Laboratório Nacional de Engenharia 
Civil, em Lisboa, tal como no caso das duas 
barragens anteriores. A Fig. 22 indica a compa- 
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ração dos deslocamentos radiais e das tensões 
médias nos paramentos dessa barragem. A Fig. 23 
refere-se à barragem de Gran Suarna e apresenta 
a comparação das tensões obtidas por modelos 
estruturais ensaiados no Laboratório Central de 
Madrid e os obtidos pelo cálculo. Trata-se de 
uma grande obra com 160 metrcs de altura. Tanto 
num caso comc noutro o acordo com as tensões 
principais máximas é excelente. Em geral, os 
valores nos pontos mais solicitados apresentam 
desvios aceitáveis, uma vez que se tome em 
consideração que os erros e as indeterminações 
da construção de modelo e das medidas das 
deformações é da ordem de grandeza de 2 + 0,02 
>< o valor da tensão, o que quer dizer que, para 
tensões da ordem de grandeza de 60 kg cm?, é 
da ordem dos 7 kg/cmº. 

Realizaram-se, últimamente, cálculos bastante 
completos da barragem de Almendra (com 200 m 
de altura), em Espanha, para várias fases de 
construção, da injecção das juntas e da pressão 
hidrostática, uma vez que a albufeira vai ser 
cheia, em várias fases, ao longo dos anos 1969 
e 1970. Os valores obtidos para os níveis 
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Fig. 21 — Barragem do SALTO DE FUNIL — Tensões principais obtidas pelo método de ajustamento 
completo e por ensaio em modelo 
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Fig. 22 — Barragem de MONTANEJOS — Deslocamentos radiais e tensões médias no eixo 
da barragem obtidas pelo método de ajustamento completo e por ensaio em modelo 


máximos concordam, muito favoravelmente, com 
resultados de ensaios de modelos conduzidos em 
1962 no Laboratório Nacional de Engenharia 
Civil. É intenção comparar esses resultados com 
observações feitas na barragem, das quais está 
encarregada a empresa CONSULPRESA, de 
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Madrid, e também com ensaios de modelo que 
estão sendo levados a cabo no Laboratório de 
Bergamo, em Itália. As Figs. 24 e 25 indicam os 
cálculos realizados e também os resultados desses 
cálculos para o nível máximo da água na albu- 
feira. 
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Fig. 23 — Barragem de GRAN SUARNA — Tensões principais obtidas pelo método de ajustamento 
completo e por ensaio em modelo 


13 — SUMÁRIO E CONCLUSÕES 


O método de ajustamento completo, para o cál- 
culo das barragens abóbada, baseia-se na divisão 
da estrutura em dois sistemas de elementos, 
um sistema de arcos e outro de consolas. Tal 
como o método «Trial load», do Bureau of Recla- 
mation, dos Estados Unidos, este método toma 
como base as hipóteses da resistência de mate- 
riais, nomeadamente de que as tensões têm uma 
distribuição linear ao longo da espessura e que 
as tensões normais radiais podem ser despreza- 
das. Em resumo, pode dizer-se que o método 
segue os princípios básicos do método «trial load», 
ainda que a resolução da distribuição de cargas 
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entre os dois sistemas de elementos se faça por 
sistemas de equações. 

A divisão da barragem em aduelas, pertencen- 
tes tanto a arcos como a consolas, permite uma 
consideração correcta das condições de equilíbrio. 
Por outro lado, consideram-se todas as forças 
interiores actuando nas faces dessas aduelas, 
assim como os deslocamentos verticais e as rota- 
ções das duas normais às faces interiores das 
aduelas. Uma vez que tanto os deslocamentos 
tangenciais como os deslocamentos verticais estão 
relacionados com as tensões de corte tangenciais, 
esses deslocamentos não podem ser considerados 
separadamente e devem-no ser em conjunto. 
As condições de continuidade e de equilíbrio 
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